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Abstract 



Static and dynamic frictional phenomena at the interface with random 
impurities are investigated in a two-chain model with incommensurate struc- 
ture. Static frictional force is caused by the impurity pinning and/or by the 
pinning due to the regular potential, which is responsible for the breaking of 
analyticity transition for impurity-free cases. It is confirmed that the static 
frictional force is always finite in the presence of impurities, in contrast to 
the impurity-free system. The nature of impurity pinning is discussed in con- 
nection with that in density waves. The kinetic frictional force of a steady 
sliding state is also investigated numerically. The relationship between the 
sliding velocity dependence of the kinetic frictional force and the strength of 

impurity potential is discussed. 
81.40.Pq, 46.30.Pa 
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I. INTRODUCTION 



Friction is an old but unsettled problem in physics [0. The Coulomb- Amonton's law of 
friction is one of the most well-known laws in physics. It has been examined that the law 
holds well under usual conditions, but breaks in some cases. It is still unsuccessfull, however, 
to clarify the conditions of Coulomb- Amonton's law to hold and the physical mechanism of 
friction. On the other hand, recent experimental measurement techniques of frictional phe- 
nomena have made great progress. One of the direction of such studies of friction tends to 
investigate nano-scale frictional phenomena. For example, atomic and frictional force micro- 
scope measurements at nanometer scale and quartz-microbalance measurements of frictional 
phenomena of monolayers have been performed extensively Stimulated by these ex- 

periments, theoretical approaches based on a microscopic model have been investigated 
extensively ^j. 

As a simple microscopic model of friction, the Frenkel-Kontorova model (FK model) 
has been employed in several works |[7|-[T0|. The FK model consists of atoms interacting with 
each other via a harmonic force under a periodic potential, which is originally sinusoidal. 
In the case that the ratio between the mean atomic spacing and the period of the potential 
is irrational, i.e., in an incommensurate case, the FK model shows a phase transition that 
is called the breaking of analyticity transition or Aubry transition 0] . When the amplitude 
of sinusoidal potential is smaller than a certain critical value, the maximum static frictional 
force vanishes. Otherwise it becomes finite. The existence of this phase transition was 
originally pointed out in a one- dimensional incommensurate FK model by Aubry [0. In 
three-dimensions, the peculiar frictional nature of the incommensurate FK model was found 
by Hirano and Shinjo They pointed out the possibility of the appearance of a super- 

lubrication state, that is the vanishing static frictional force state, between incommensurate 
surfaces in real materials. Matsukawa and Fukuyama investigated the static and kinetic 
properties of friction in a one-dimensional model |jTl|. The model employed in their study 
consists of two atomic chains, where interchain force works between atoms in one chain and 
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those in another and harmonic force works between neighber atoms in each chain. The effect 
of energy dissipation and the relaxation of atoms in both chains are taken into accout prop- 
erly. On the basis of their model, a definitive expression of both static and kinetic frictional 
force can be derived. They found that the critical value of the potential strength of the 
breaking of analyticity transition depends strongly on the rigidity of the chains. In the case 
that atoms in the lower chain are fixed periodically, which corresponds to the FK model, 
that critical value is large. When atoms in the lower chain can relax, that critical value 
decreases and the parameter range for the state with vanishing maximum static frictional 
force becomes small. They also found that the sliding velocity dependence of the kinetic 
frictional force becomes weaker as the maximum static frictional force increases. 

In these theoretical studies of friction, surface randomness has not been taken into ac- 
count. In a realistic system, however, there always exists surface randomness. The effects of 
surface randomness should be taken into consideration to understand frictional phenomena 
in realistic systems. Based on a variational method and a perturbational analysis, Mat- 
sukawa and Fukuyama showed that the state with vanishing maximum static frictional force 
becomes unstable against surface randomness The frictional transition is smeared out 
because the randomness breaks the translational symmetry of the system. Hence a finite 
maximum static frictional force appears even if the interatomic interaction is sufficiently 



small [|T3 |. 



While finite static frictional force is expected to appear in systems with surface ran- 
domness, the magnitude of the maximum static frictional force has not been clarified yet. 
Furthermore, surface randomness influences the kinetic friction. The velocity dependence of 
the kinetic frictional force would be different from that of the system with clean surfaces. 
In this paper, using a one-dimensional incommensurate model, we investigate the static and 
kinetic frictional forces subject to surface randomness. 

This paper is organized as follows. In section II we define a one-dimensional incom- 
mensurate model of friction and introduce surface randomness. Based on a numerical study 
and the impurity pinning theory of density waves, we discuss the static frictional force in 



section III. The kinetic frictional force, in particular its velocity dependence, is investigated 
in section IV. Section V is devoted to summary and discussion. 



II. THEORETICAL MODEL OF FRICTION 

In order to calculate frictional force, we consider a one-dimensional model of friction 
proposed by Matsukawa and Fukuyama . This model can deal with the friction between 



two deformable atomic chains interacting each other. To simplify the problem, however, 
we here assume that one of the two chains is fixed and the other chain can be deformed. 
The atoms of the fixed chain are periodically arranged at a regular position. We call the 
deformable chain the upper body and the fixed chain the lower body. The atoms in the 
upper body have a one- dimensional degree of freedom in the direction parallel to the chain. 
The interatomic force of the upper body is approximated by the harmonic one. The effect 
of the energy dissipation is taken into account by the term in the equation of motion which 
is proportional to the difference between the velocity of each atom and that of the center 
of gravity. The external force is applied to the upper body to the horizontal direction. 
Assuming an overdamped condition, we get the equation of motion of the i-th atom of the 
upper body as follows. 

rrialaiUi - {Ui)i} 

= Ka{Ui+i + Ui-i - 2Ui} + Fi{Ui - Vj) + Fex, (l) 

jeb 

where Ui (fj) is the position of the i-th atom of the upper (lower) body, m the atomic 
mass, 7 the parameter of energy dissipation, Ka the strength of the interatomic froce, 
and {ui)i represents the average of iii with respect to i. Fj and F^x are the interchain 
force between atoms in different chains and the external force, respectively. We adopt the 
following interchain potential: 



Ui = — — exp 



(2) 
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where Ki is the interchain potential strength and the mean atomic spacing of the lower 
body. The interchain force is given by Fi{x) = —-^Uj. Since we now consider the case 
where all the atoms of the lower body are fixed, the total interchain force that acts on the 
i-th atom of the upper body is a function of Ui only, which is obtained by summing up the 
contribution from all the atoms of the lower body. 

J:Fi(u,-Vj)^F,(u,). (3) 

This means that the present model corresponds to the kinetic FK model with energy dissi- 
pation. 

Now wc consider surface randomness such as impurities. It is assumed that the interchain 
potential consists of two parts, i.e., the original regular potential and the impurity potential, 
which is distributed randomly. Furthermore we assume that the impurity potential has the 
same form as the original one except for the potential strength. Hence the total interchain 
potential is given by 

Nt Ni 

Ui^{ui) = ^Ui{ui-vj)+ Yl Ur^{ui-vj), (4) 

j=l j=impurity sites 

where U[ represents the original interatomic potential, C//™^ the impurity potential and 
is the number of impurities. The impurity potential, t//™'^, is written as 

where Kf"^'^ is the strength of the impurity potential. In this model, the frictional force 
Ffric in a steady state can be calculated by the following equation: Ff^c = — {Fj{ui))^ , 
where (. . .)^ means the temporal average and Fj{x) — —-^Ui^{x). Because the system is in 
a steady state Ffric is equal to the total external force, NaFex- 

The equation of motion is numerically solved using the Runge-Kutta formula. In the 
calculation the periodic boundary condition, Ui — Ui+N + L, is adopted, where L is the 
system size. Then the ratio Ca/ci, is equal to the ratio N^/Na, where Ca is the mean lattice 
spacing of the upper body and Na and are the number of atoms of the upper and lower 
bodies, respectively. 
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x) = 



exp 



We here consider an incommensurate system with periodic boundary condition and de- 
termine the ratio Ca/cb based on the continued-fraction expansion of the golden mean. 

Ca Nb 144 , , 

— = -f = = 1.6179- ■■, 6 

where the expansion is truncated at the tenth order. The parameters of the present model 
are set as 

144 

Ca = —,Cb = l,m = Ka = ^ = l,Nj = 70. (7) 
oy 

The number of impurities Nj is fixed throughout the present numerical simulation. With 
increasing external force, the maximum static frictional force, i.e., the threshold force at 
which the upper body starts to slide, is numerically estimated. The kinetic frictional force 



is evaluated in steady sliding states with constant velocities [11 



III. MAXIMUM STATIC FRICTIONAL FORCE AND IMPURITY PINNING 

EFFECTS 

For convenience of the following discussion, we define some energy scales that characterize 
the present model. The impurity potential energy, the elastic energy, and the coupling energy 
to the external force per impurity are expressed as Ei, Eei, and Ef, respectively, and are 
given by 

E, = ^Xr^, (8) 

Eel = \Ka X Ni/Nb, (9) 

Ej = Fe, X Nb/Ni. (10) 

For the comparison between these energy scales it is convenient to introduce the parameters 
Si and Sf as follows: 

Si = EjEei, (11) 
Sf = Ef/Eei. (12) 
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The strengths of the impurity potential and the external force are discussed in terms of Ei 



and Ef, respectively. According to the impurity pinning theory of density waves [p^|-p!6[] the 
nature of the impurity pinning is quite different whether Si ^ 1 or Si <^ 1. 

We now attempt to interpret the impurity pinning effect as the origin of static friction. 
The scenario is as follows. For large the effect of the impurity potential overcomes those 
of the elastic energy, and strong pinning states, in which each impurity pins the chain and 
the large lattice distortion is accompanied, are preferred. Maximum static frictional force is 
expected to be proportional to the pinning energy, which is the sum of the impurity potential 
energy and the elastic energy in the pinned state. Hence, the normalized maximum static 
frictional force ej^'^' behaves as 

ef"- oc Ei. (13) 

On the other hand, for small values of Ei each impurity cannot pin the upper body due to 
the resulting large loss of elastic energy. Instead collection of impurities in a finite domain 
pin the upper body. The size of the domain is determined so as to minimize the pinning 
energy. In each domain the system can gain the impurity energy in the scale of fluctuation. 
Then we get 

ef^'-ocEA (14) 

in the case of one-dimensions. 

The consideration above, however, does not take into accout the effect of the originally 
underlying regular potential Ui. Even for impurity- free cases there can exist the pinning 
behavior in the present system, i.e., the appearance of the finite maximum static frictional 
force due to the Aubry transition for Kj larger than a certain critical value. This effect is 
crucial in considering the static frictional force. 

Figs. l(a)-(c) show eJ-'^' versus Ei curves for the several strengths of the interchain 
potential Kj, where the maximum static frictional forces obtained are averaged over 3 to 
10 samples. For Kj = 0.05 (Fig. 1(a)), the weak impurity pinning relation (eJ'^' oc Ei^) 
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and the strong impurity pinning relation (s™-^- oc Si) hold well in the regimes Ei < 10 and 
£j > 20, respectively. There exists a crossover point between the two characteristic impurity 
pinning states around Si ^ 20. Additional crossover behavior appears for Kj = 1 where the 
corresponing clean systems have finite static frictional force. As observed in Fig. 1(b), the 
crossover between weak and strong impurity pinning states found for Kj — 0.05 apparently 
exist around Si ~ 20 even for Kj = 1. The regime of the weak impurity pinning states 
for Ki — 1, however, is rather narrowed, and apparent deviation from the weak impurity 
pinning behavior is observed for < 2. This deviational behavior is considered a crossover 
between the weak impurity pinning state and an another state. The latter is the state pinned 
mainly by the regular potential and is close to the breaking of analyticity state due to the 
Aubry phase transition in the impurity-free case. In this regime the impurity potential is 
considered approximately to be a weak perturbation to the impurity-free state. However, 
in the regimes > 2, the effect of the impurity potential becomes dominant, and then the 
weak pinning state appears. The deviation point from the weak pinning state shifts to the 
large £j regime as Kj increases. For a strong interaction {Kj = 6, Fig. 1(c) ), the weak 
pinning regime is not observed. It is confirmed from the fact that the exponent of eJ^'^^ on 
£j in the regime £j < 10 is quite different from that of the weak pinning state. The weak 
pinning states are destroyed by the regular interchain potential because this potential always 
dominates the impurity potential in the weak impurity pinning regimes £j < 10. 



IV. KINETIC FRICTIONAL FORCE 



In this section we investigate the effect of the impurity potential on the velocity de- 
pendence of the kinetic frictional force. Figs.2(a-b) show the velocity(i') dependence of the 
kinetic frictional force for various strengths of the impurity potential. These values of the im- 
purity potential correspond to the weak pinning regime. In addition to the finite Si cases we 
show the kinetic frictional force of the clean system(£i = 0) for comparison. In a weak regular 
interchain interaction case, Kj — 0.05(Fig.2(a)), it is obviously observed that the kinetic fric- 



8 



tional force for Ei = exhibits velocity strengthening(Ffric oc v) and weakening (Ffric oc v~^) 
features in low and high velocity regimes, respectively. Such behavior is excellently explained 
by a perturbational analysis In the case of finite Ei, the behavior is quite different 

and the frictional force tends to a finite value in the limit of vanishing velocity [0. It is an 
effect of the impurity potential. For small the effect is more significant at low velocities 
than at high velocities, and additional velocity- weakening kinetic frictional force appears in 
the low volocity regime. As the strength of the impurity potential increases the maximum 
static frictional force becomes larger than the kinetic frictional force. We here note that 
the behavior of kinetic frictional force in disordered systems was discussed theoretically by 



Sokoloff pO|. When applying his method to one-dimensions, a velocity- weakening feature 
of kinetic frictional force can be deduced based on the perturbation theory for a general 
disordered potential. The result is, however, obtained in the underdamped regime, while 
the present results are obtained in the overdamped regime. Fig. 2(b) shows the result on the 
kinetic frictional force for Kj = 1 where there exists large maximum static frictional force. 
In this case the kinetic frictional force for Ei = shows a similar velocity dependence to that 
for finite Ei, that is, almost velocity-independent behavior in the low- velocity regime and 
inversely proportional behavior to the sliding velocity in the high-velocity regime. In both 
cases of Figs. 2 (a) and (b) the strength of the kinetic frictional force monotonically increases 
with increasing Ei. 

The kinetic frictional force in the sliding states evolved from strong impurity pinning 
states is shown in Fig.3. Under this condition, a very large maximum static frictional force 
appears and the static and kinetic friction is mainly caused by the interaction between the 
chain and the impurities. The kinetic frictional force in Fig.3 is similar to that in Fig.2(b), 
i.e., the velocity-independent and velocity- weakening features. No distinct qualitative differ- 
ence in the behavior of the kinetic frictional force is observed between the different strengths 
of the interchain potential, Kj = 0.05 and 1. 

From the above results on kinetic frictional force, we found that the velocity-weakening 
behavior in the high-velocity regime is a common feature of the present model and the 
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velocity-independent behavior in the low-velocity regime is peculiar to the system with 
large maximum static frictional force. 



V. SUMMARY AND DISCUSSION 

We have calculated the static frictional force by computer simulations using a one- 
dimensional incommensurate model of friction. It has been numerically confirmed that the 
randomness due to impurities destroys the vanishing static frictional force states and the 
maximum static frictional force always becomes finite. This is consistent with the variational 
calculation by Matsukawa and Fukuyama 0]. 

We have also found that frictional phenomena have a close relationship to impurity pin- 
ning phenomena in density wave systems. For the case of weak regular interchain interaction 
the nature of impurity pinning states changes from weak pinning to strong pinning as the 
strength of impurity potential increases. Similar behavior is also observed in density waves 



I^ , |2T[] . For a strong regular interchain interaction the maximum static frictional force is 
finite even in impurity-free systems. The feature of the impurity pinning in the weak pin- 
ning regimes is modified. In this case the pinning is caused essentially by the strong regular 
potential rather than by the impurity potential. The crossover behavior of the maximum 
static frictional force has been found from the pinned sate affected by the regular potential 
to the weak or strong impurity pinning state as the impurity potential is strengthened. 

We have also investigated the kinetic frictional force in the presence of impurities. It 
has been found that surface randomness is crucial to the velocity dependence of the kinetic 
frictional force. The behavior of the kinetic frictional force in the present work should be 
compared with that of the Coulomb- Amonton's law, i.e., the kinetic frictional force is less 
than the maximum static frictional force and does not depend on the sliding velocity. In the 
case of large maximum static frictional force the velocity-independent behavior of kinetic 
frictional force is observed in the low-velocity regime. The magnitude of the kinetic frictional 
force is, however, equal to the maximum static frictional force. 



10 



In a realistic system there are several complicated factors that are not taken into consid- 
eration in the present theoretical model. For example, the effect of other types of randomness 
such as steps on surfaces is not examined here. Such an effect may be also an important 
factor, but we expect that the velocity dependence of the kinetic frictional force would not 
depend on the details of the randomness and would show a similar behavior to that of the 
present model because the major contribution to the kinetic frictional force comes from the 
excitation of phonons with long wavelength. Nevertheless several important issues are left 
unsolved and therefore we need further investigation to compare the detailed features of the 
friction model with those of realistic systems. 
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APPENDIX: 

We make some comments on Eqs. (13) and (14). As mentioned in the text, these 
relations are derived using the impurity pinning theory of density waves. Concerning the 
analogy between the present model and the model of density waves, the following points 
should be noted. In the case of density waves [0, the average impurity energy vanishes 
when density waves are rigid and have perfect periodicity. The energy gain due to the 
impurity potential arises from the distortion of density waves. On the other hand, since the 
present model has the attractive potential t/}™^, the energy gain always exists even without 
the distortion of the chain. This energy, however, does not affect the frictional force because 
the amount of the energy is independent of the position of the chain. In both systems the 
energy cost caused by the distortion is relevant to the pinning effect. Hence, in the present 
model, the same behavior of the impurity pinning as that in density waves is observed. 
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Fig. 1. Si versus ef'^\ The solid and dashed lines represent the relations: ej'^' oc £j3 and 
ef-^- oc Ei, respectively. When Ei = 0(Xi"P = 0), ef^- = 0, 0.760, and 18.0 for Kj = 0.05, 1, 
and 6 respectively. 

Fig. 2. The velocity dependence of kinetic frictional force for weak pinning, (a) 
Kj — 0.05. Squares, circles, diamonds, triangles, and crosses represent the results for 
Ei = 0.926, 0.309, 0.0617, 0.0206, and {K'^^ = 0.45, 0.15, 0.03, 0.01, and 0), respectively (b) 
Ki — 1. Squares, circles, triangles and crosses represent the results for £j = 18.5, 6.17, 0.411, 
and (X™^ = 9, 3, 0.2, and 0), respectively. Arrows indicate the maximum static frictional 
force. 
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Fig. 3. The velocity dependence of kinetic frictional force for strong pinning. Circles and 
triangles represent the results for Si = 51.3 {K\"''^ = 24.95) and Kj — 0.05 and for £j = 90.5 
j-^imp _ ^^-j _ -^^ respectively. Arrows indicate the maximum static frictional force. 
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Fig. 1(c) Kawaguchi and Matsukawa 
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